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In this paper we analyze nuclear medium effects on DVCS amplitudes in the Xbj range of 10~^ — 
10~* for a large range of and four different nuclei. We use our nucleon GPD model capable of 
describing all currently available DVCS data on the proton and extend it to the nuclear case using 
two competing parameterizations of nuclear effects. The two parameterizations, though giving 
different absolute numbers, yield the same type and magnitude of effects for the imaginary and real 
part of the nuclear DVCS amplitude. The imaginary part shows stronger nuclear shadowing effects 
compared to the inclusive case i.e. ^2^, whereas in the real part nuclear shadowing at small xi,j and 
anti-shadowing at large Xbj combine through evolution to yield an even greater suppression than 
in the imaginary part up to large values of xi,j. This is the first time that such a combination of 
nuclear effects has been observed in a hadronic amplitude. The experimental implications will be 
discussed in a subsequent publication. 



PACS numbers: ll.lO.Hi, ll.30.Ly, 12.38.Bx 



I. INTRODUCTION 

Currently information about parton structure of nu- 
clei is very limited. Most of the data were obtained at 
rather large Xbj > 5 x 10^^ where deviations of the nu- 
clear structure from additivity for F2A are of the order of 
a few percent, until one reaches the region of the EMC 
effect Xhj > 0.4 were the structure functions themselves 
are already pretty small. For gluons the situation is even 
more uncertain - the theoretical analysis based on QCD 
momentum sum rules and the expectation of a significant 
nuclear shadowing P, Q leads one to expect a gluon en- 
hancement of the order of 10 — 20% for Xbj ~ 0.1. The 
analysis of the NMC data supports such a conclusion 
3, 4|. However, the range of the NMC data is rather 
limited and higher twist effects may affect the analysis 
substantially. 

The interpretation of the data at smaller Xbj where 
one observes the shadowing of F2A, is even more difficult. 
This is due to the strong correlation of Xbj and in the 
relevant kinematic region of x^j 

In this kinematic region theoretical calculations of shad- 
owing based on the Gribov theory connecting diffraction 
in eN scattering and shadowing in eA scattering describe 
the data well. However, in these analysis the vector me- 
son contribution which is a higher twist effect was found 
to be substantial thus indicating the importance of higher 
twist shadowing effects [ElEH- principle, local du- 
ality could result in the vector meson contribution being 
dual to the continuum. However, the recent analysis in 8] 
based on the theory of leading twist nuclear parton shad- 
owing P, 1^ has explicitly demonstrated that 30 - 40% 
of the shadowing for F2A{xbj < 10~^,Q'^ < 2 GeV^) is 
due to higher twist effects. 

Clearly in this situation one needs new experimen- 



tal tools to study the nuclear partonic structure. Ex- 
clusive reactions appear to be a very natural candidate. 
Among them deeply virtual Compton scattering (DVCS) 
(see Fig. is the simplest reaction since the final state 
real photon cannot be involved in the rescatterings, while 
in the case of meson production a suppression of final 
state interactions (color transparency) probably requires 
rather large values of Q^. 

The most interesting production channel to investi- 
gate is the one where the nucleus remains intact - co- 
herent nuclear DVCS. The generalized parton distribu- 
tions (GPDs) which enter in this case carry much more 
direct information about the nuclear structure than the 
channels with a nuclear break-up. Besides, the extrac- 
tion of the coherent signal is much easier experimentally 
both for fixed target experiments with relatively modest 
resolution on the missing mass, and for future collider ex- 
periments like the planned Electron Ion Collider (EIC) 
where it will be possible to use a zero angle calorimeter 
to select events were the nucleus remained intact. 

This, however, comes at a price. The average values 
of t in the coherent scattering are small ~ and can 

hardly be measured. Hence, one has to deal with cross 
sections integrated over t. The challenge in this case will 
be to observe the DVCS signal at small t and small Xbj 
where the QED Compton or Bethe-Heitler process is very 
important. 

Thus our task in studying nuclear DVCS is two fold. 
First, we want to construct generalized nuclear parton 
distribution functions relevant for the discussed process 
and then to calculate the DVCS amplitudes beyond the 
impulse approximation which was considered in Til for 
both coherent and incoherent channels. This will be the 
subject of the first of two papers. In the second paper we 
will study various DVCS observables to find out which of 
them are feasible for small Xbj DVCS studies. 

Our strategy will be to combine techniques which 
we developed in the modeling of next-to-leading order 
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(NLO) GPDs for nucleons in Il3 and which describe the 
current world DVCS data [13 with fairly high accuracy 
and the recent calculations of nuclear PDFs at small Xbj 
iniM^I' These calculations are based on the connection 
between diffractive parton densities measured at HERA 
and shadowing for nuclear parton densities and indi- 
cate that at small Xhj shadowing should be large in both 
the quark and gluon channels. 

There are several aspects of nuclear DVCS at small 
Xbj which appear interesting. One goal would be to 
test our understanding of the connection between PDFs 
and GPDs. Qualitatively, we expect that the shadow- 
ing regime will start for the imaginary part of DVCS 
at larger Xbj since the effective Xbj which enters in the 
DVCS amplitudes is of the order of Xbjf^. The space- 
time picture of DVCS points in the same direction since 
the coherence length in the transition of an intermediate 
state of mass, M, into a real photon is larger than in the 
case of DIS - the non-covariant propagator is ^ 1/M^ 
instead of 1/{M'^ + Q^)- Another aspect is the impor- 
tance of the gluon GPD feeding into the quark singlet 
GPD in the course of evolution especially in the NLO 
approximation which was demonstrated in Since 
the shadowing of gluons and quarks turns out to be dif- 
ferent within the model adopted in this paper, we expect 
the A-dependence of the DVCS observables to be differ- 
ent from that of F2A ■ For a comparison we will give also 
results based on the model of which is based on a fit to 
Xbj > 10~^ data assuming leading twist dominance and 
extrapolating to smaller Xbj with the assumption that 
the gluon and quark parton densities are screened in the 
same way at the initial scale of perturbative QCD evo- 
lution. Another unique aspect of DVCS is the ability to 
measure the real part of the amplitude. In this case we 
expect significant nuclear effects already at rather large 
Xbj ~ 0.1 since the real part of the DVCS amplitude is 
sensitive to the rate of variation of parton densities with 
Xbj, which, for the nuclear case, is already rather different 
at Xbj ^ 5 X 10^^. 

Studies of nuclear DVCS will have important spin offs 
for studies of color transparency in coherent high energy 
reactions since the same GPDs enter in both DVCS and 
vector meson production. Hence comparing results for 
these two processes will allow one to find the value of 
at which perturbative color transparency becomes essen- 
tially model independent. 

This paper is organized as follows: In section ^] we 
summarize the procedure for modeling nuclear GPDs 
based on the information about nucleon GPDs. We will 
also discuss the accuracy of factorizing GPDs as a prod- 
uct of the GPD at i = and the nuclear form factor. 

In the next section we summarize the previous findings 
on NLO DVCS amplitudes and give numerical results for 
a set of nuclei. We will demonstrate that nuclear effects 
are indeed quite strong and different from the case of F2a- 
We find that there is a significant difference between nu- 
clear effects between the two cases. Interestingly enough 
we find that the nuclear effects for the real part of the am- 




FIG. 1: (Color online) a) DVCS graph, b) BH with photon 
from final state lepton and c) with photon from initial state 
lepton. 



plitude may extend all the way to the highest Xbj where 
the coherent peak is still present, that is Xbj ~ 0.2. This 
is due to a curious interplay between nuclear shadowing 
and anti-shadowing for 10^^ < Xbj 

In the second paper we will estimate various DVCS 
observables. The main focus will be on the kinematic re- 
gion of the EIC collider, though we will also give numbers 
for fixed target kinematics. Overall we find that DVCS 
signals are large enough for all standard observables to 
be measured with high accuracy at the EIC. 



II. A NUCLEAR GPD MODEL AND NUCLEAR 
DVCS AMPLITUDES 

In the following, we will use and review the model 
for twist-2 GPDs first introduced in 12] using the off- 
diagonal representation J^J of GPDs. 

Based on the aligned jet model (AJM) (see for example 
[Tell) the key Ansatz of [l3| in the DGLAP region is: 
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(1) 



where refers to any forward distribution and 
HS-NS^XX) = Hi{X,0 ± H<i{X,C)- This Ansatz in 
the DGLAP region corresponds to a double distribution 
model ^17', Ts, Is] with an extremal profile function al- 
lowing no additional skewdness save for the kinematic 
one. 

Note that we choose a GPD representation first intro- 
duced in |l5l |. which is maximally close to the inclusive 
case i.e A e [0, 1], C = Xbj with the partonic or DGLAP 
region in [C, 1] and the distributional amplitude or ERBL 
region in [0, C]- Note here that in the case of nuclei we 
speak of the momentum fractions per nucleon, therefore, 

Ce [0,A]. 
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The prescription in Eq. does not dictate what to 
do in the ERBL region, which does not have a forward 
analog. The GPDs have to be continuous through the 
point X = and should have the correct symmetries 
around the midpoint of the ERBL region. They are also 
required to satisfy the requirements of polynomiality: 
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with X = ■ The ERBL region is therefore modeled 

with these natural features in mind. One demands that 
the resultant GPDs reproduce the first moment Mi = 3 
and the second moment M2 = 1 + CC/{2 - C)^ Jifl. C 
was computed in the chiral-quark-soliton model |2ll| and 
found to be —3.2 and is related to the D-term |23| which 
lives exclusively in the ERBL region. This reasoning sug- 
gests the following simple analytical form for the ERBL 
region {X < (): 

H9'^S{X, = iJf'^^(C) [1 + As^^^iQC'^'^'^iX, 0] , 



H^iX,0 = H'{C) 



C/2 



[1 + A^(C)C^(X,C)] , 
(3) 



where the functions 
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vanish at X = C to guarantee continuity of the GPDs. 
The A*(C) are then calculated for each C by demanding 
that the first two moments of the GPDs are explicitly sat- 
isfied. For the second moment, what is done in practice is 
to set the D-term to zero and demand that for each fiavor 
the whole integral over the GPD is equal to the whole in- 
tegral over the forward input PDF without the shift. For 
the final GPD, of course, the D-term is added to the 
quark-singlet (there is no D-term in the non-singlet sec- 
tor) using the results from the chiral-quark-soliton model 
|2ll |. The gluonic D-term, about which nothing is known 
save its symmetry, is set to zero for Qq- Due to the 
gluon-quark mixing in the singlet channel, there will be 
a gluonic D-term generated through evolution. 

It would be straightforward to extend this algorithm 
to satisfy polynomiality to arbitrary accuracy by writing 
the A*(C) explicitly as a polynomial in C where the first 



few coefiicients are set by the first two moments and the 
other coefficients are then either determined by the ar- 
bitrary functional form, as is done here, or, perhaps the- 
oretically more appealing, one chooses orthogonal poly- 
nomials, such as Gegenbauer polynomials, for which one 
can set the unknown higher moments equal zero. Phe- 
nomenologically speaking, the difference between the two 
choices is negligible. 

The above Ansatz also satisfies the required positivity 
conditions 0,113,13 and is in general extremely flexible 
both in its implementation and adaption to either other 
forward PDFs or other functional forms in the ERBL 
region. Therefore, it can be easily incorporated into a 
fitting procedure making it phenomenologically very use- 
ful. Since we will concern ourselves only with NLO, in 
what follows we will use CTEQ6M 25] as the forward 
distribution. 

In order to extend this parameterization for the proton 
to a nucleus, we have to include the effects due to the 
nuclear medium. We do this by multiplying the GPD in 
Eq. ^ with the ratio of nuclear per nucleon to nucleon 
inclusive parton distribution function (PDF) shifted by 
the same amount as the forward distribution in Eq. 
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The ratio of nuclear to nucleon PDF is taken from 
. Note that the parametrization of Eq|Sl corresponds 
to the same value of Mi as for the nucleon (which is any- 
way necessary to satisfy the baryon charge conservation) 
and also leads to an A-independent D-term. Note that 
it is difficult to generate an A-dependent D-term via an 
attachment of the current to the simplest t-channel ex- 
changes. As an example consider the contribution to the 
D-term due to two pions or a resonance in the t-channel. 
The two-pion term where the pions are attached to two 
different nucleons, is strongly suppressed for large A since 
it requires a spin fiip of both nucleons. Furthermore, it 
is quite difficult to attach a resonance to two nucleons. 
At the same time, one cannot exclude the possibility of 
a coherent interaction of a current with several nucleons. 
In fact, the real part of the amplitude appears to be more 
sensitive to large distance effects which may enhance co- 
herent effects. Hence, one cannot exclude that the D- 
term actually does depend on the atomic number [soj. 
However, the quantities we consider are not very sensi- 
tive to the D-term and so taking the D-term to be inde- 
pendent of A seems to be a natural first approximation. 

As far as the f-dependence is concerned, we postpone 
the discussion to the next paper dealing with nuclear 
DVCS observables. Thus we will consider the case t — 
for the time being. This is sufficient at present, since 
larger values oi t » tmm will not be experimentally 
accessible in the coherent peak we are interested in. 

After having evolved the GPDs using the same pro- 
gram successfully employed in the real and imagi- 
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nary part of the twist-2 DVCS amplitude in LO and NLO 
given below are calculated using the same program as in 
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where stands for the nuclear GPD Ha and for the 
respective Wilson coefficient functions. 

The +ie prescription is implemented using the Cauchy 
principal value prescription ("P.V.") through the follow- 
ing algorithm: 

p.v. dxri^- 1^ F{x, c, Q') = 

j^^ dXT(^-^- 1^ {T[X, c, - ^(C, C, Q')) + 
+ ^(C,C,Q')^'rf^T(^-i) • (8) 

The relevant LO and NLO coefficient functions are the 
same as for the nucleon and can be found in l2a, l23 ■ 



III. 



RESULTS FOR NUCLEAR DVCS 
AMPLITUDES 



In the following we give our results for the unpolarized 
nuclear DVCS amplitudes in NLO for the following nu- 
clei: O- 16, CA~40, Pd~ 110 and P6-208. We will use 
the FGS sj and the Eskola et al. 3 parameterizations 
of nuclear medium effects. The initial scale for the per- 
turbative evolution in both cases will be Qg = 2 GeV^. 
We choose such a low scale to be able to make state- 
ments about HERMES kinematics (0.05 < Xbj < 0.3 and 
1 < Q2 < 9 GeV^) in our forthcoming paper on nuclear 
DVCS observables. The FGS parameterization has an 
initial scale of Qq = 4 GeV^ and we, therefore, evolve it 
backward to our Qq. The parameterization of Eskola et 
al. has a Qq = 2.25 GeV^ and is thus close enough to our 



starting scale. Unfortunately, it is a LO parameterization 
which we will use in NLO. 

Our treatments of these two parameterizations can be 
rightfully called into question. However, our aim at this 
point is not to make precision statements about nuclear 
DVCS but rather look for the type of effects which can be 
expected and how parameterization dependent they are. 
If the differences in parameterizations will also material- 
ize in DVCS observables in a noticeable way, experiments 
we will have a good chance of discriminating between dif- 
ferent parameterization relying on different assumptions 
of the underlying physics. 

Let us first concentrate on the imaginary part of 
the unpolarized nuclear DVCS amplitude which depends 
only on the DGLAP region. In Fig. [3 we show the ratio 
of the imaginary part of nuclear to nucleon DVCS ampli- 
tude in NLO for fixed and varying Xbj and in Fig. |2| 
we show the same ratio for fixed Xbj and varying Q^. 

As is easily seen from Fig. |21 the nuclear shadowing 
corrections differ substantially between the two parame- 
terizations at low and are very similar for larger Q^. 
However, the nuclear shadowing corrections are in both 
cases stronger than in the forward case of F.f/F2 (see 
for example 8]). In particular in the region of small 
Xbj < 5 X 10~^ one observes a significant slope increase 
in Xbj compared to the forward case. This behavior is 
easily explained with first the shifted input toward lower 
values of X. This leads to enhanced nuclear shadowing 
in the small Xbj region since shadowing continues to in- 
crease toward smaller values of Xbj- Secondly, since the 
value of the imaginary part of the amplitude is mainly 
determined by the region X ~ 0{Q of large light-like 
distances (see and since this region is enhanced by 
perturbative evolution, one expects that until the nuclear 
enhancement effects from large X are degraded enough 
in momentum through evolution, nuclear shadowing ef- 
fects will increase with the increase in Q^. As can be 
seen from Fig. O evolution for fixed Xbj will indeed 
initially increase nuclear shadowing but the vanishing 
of nuclear shadowing with an increase in as seen in 
/ F2 does not happen as rapidly as in DVCS. This oc- 
curs precisely because the relevant region for the value 
of the amplitude is the one of large light-like correlations 
X ~ 0(C) which is enhanced through perturbative evo- 
lution. Large-light like correlations in turn are more sen- 
sitive to nuclear medium effects and thus nuclear shad- 
owing will remain stronger for a wider range in . We 
also observe a stronger A dependence at small than at 
large Xbj which is again stronger than the forward case. 
Once more the shift to relative smaller values of Xbj is the 
reason. Nuclear shadowing is stronger for Ph — 206 than 
for — 16 and thus a relative shift to smaller values of 
Xbj will increase this difference further. Evolution does 
not change this picture a lot since the evolution effects 
are the same for all nuclei. As we just explained, nuclear 
shadowing persists under evolution for a wider range in 
than in the forward case where nuclear shadowing 
disappears more rapidly and hence differences between 



5 



R, =ImA 7ImA,, in NLO R, =ImA7ImA,, in NLO 

Im A N Im A N 




0.6- 



0.4- 



0.2- 



Q' = 4Ger 
I I I — I I I — rm 



— 0-16 
■■■■ Ca-40 
Pd-110 
■-■ Pb-208 



qI I I I I I I 




0.6- 



0.4- 



0.2- 



Q" = 10GeV" 

I I I 1 I I TTTT 



R, =Im7Im„inNL0 

Im A N, 
X =10 



x,, = 10" 



— 0-16 

■■■■ Ca-40 

-- Pd-110 

■- Pb-208 



0.7 7, 
0.6 -" 



0.5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i' 



— 0-16 
■■■■ Ca-40 
Pd-110 
■-■ Pb-208 



1 n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ij 1 

0.9 - - 0.9 



10 20 30 40 50 60 70 80 90 
QlGeV"] 
x,, = 10"" 



4 R, 

Im 

0.7 

-I 0.6 

0.5 



1 1 1 M 1 1 1 M 1 1 M 1 M 1 1 


1 M 1 M 1 1 1 M 1 1 M 1 1 1 1 


!^ 








'i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ll' 



10 20 30 40 50 60 70 80 90 
QlGeV'] 
X =0.1 



0.9 r: 



R, 



0.7 - 
0.6- 



10"^ 10"^ 10" 10' 



R, =ImA7ImAJnNL0 

Im A N 
Q' = 4GeV 



0.9- 



0.7 
0.6 h 

0.4 |- 

0.3 7 

0.2 '- 

0.1 - 








— 0-16 

■■■■ Ca-40 

-- Pd-110 

■-■ Pb-208 



l_Lli 



qI I I I I I I 

-4 -1 -1 -1 

10 10" 10" 10 



R, =ImA7ImAJnNL0 

Im A N 
Q" = 10GeV" 

I I I 1 I I TTTT 



0.5 



Uj 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I I 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I l' 



0.9 - 



'■ R 



10 20 30 40 50 60 70 80 90 
0"[GeV"] 

R. =Im./Im„inNLO 



Im 

0.7 

-I 0.6 h 

0.5 



'<r<'|'i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i' 



10 20 30 40 50 60 70 80 90 
QlGeV"] 



im A N, 



0.9- 



0.7 
0.6 h 

0.4|- 
0.3 1- 
Q.l'- 
0.1- 








— 0-16 

■■■■ Ca-40 

-- Pd-110 

■- Pb-208 



I I l_LU 



— 0-16 

■■■■ Ca-40 

-- Pd-110 

■-■ Pb-208 



0.9 - 
0.8 f 

. 

0.7 r 
0.6- 



0.5 



\ = 10 X,, = 10 

bj bj 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I I I I I I 1 1 1 1 1 1 1 1 1 1 1 1 IJ 1 Ml 1 1 1 I 1 1 1 1 I I 1 1 1 I I 1 1 1 I I 1 1 1 I I 1 1 1 I I 1 1 1 1 I l_ 

0.9 



^.1... I. ..I. ..I. ..I. ..I. ..I. ..I. ..I. 



■ R ''''' 

: Im ;• 

0.7^ 
0.6 - 



10 20 30 40 50 60 70 80 90 
Q'[GeV'l 
K = 10"^ 



0.5 



ilii I III! I I III I I ill I I III I I III I mini' 



-4 -3 -2 -1 -4 -3 -2 -1 

10 10 10 10 10 10 10 10 




10 20 30 40 50 60 70 80 90 
Q"[GeV"] 
x, = 0.1 



i 


1 1 M 1 1 1 M 1 1 M 1 1 1 l_ 


'i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



10 20 30 40 50 60 70 80 90 
QlGeV"] 



10 20 30 40 50 60 70 80 90 
QlGeV'] 



FIG. 2; (Color online) Ratio of nuclear per nucleon to nucleon 
DVCS amplitudes vs. xtj for two values of for the FGS 
(upper plot) and Eskola (lower plot) parameterizations. 



FIG. 3: (Color online) Ratio of nuclear per nucleon to nucleon 
DVCS amplitudes vs. for four values of Xbj for the FGS 
(upper plot) and Eskola (lower plot) parameterizations. 
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nuclei as well. When comparing the two parameteriza- 
tions one notices the peculiar behavior of the ratio for 
the FGS parameterization at — A GcV^ in the region 
of 5 X 10"'^ < Xbj < 10"^. This is due to the fact that in 
the FGS parameterization the ratio of nuclear to nuclcon 
ghion PDF is above 1 in the region 10^^ < Xbj < 10^^ 
and we are using shifted distributions as our input for the 
GPDs. As it turns out the imaginary part of gluon ampli- 
tude is negative whereas it is positive for the quarks |2g 
and the relative suppressions in the quark and the gluon 
are sufficiently different enough in nuclei and nucleons 
that in the sum this difference disappears! Note that this 
is not the case for the Eskola parameterization. However, 
at = 10 GeV^ this behavior has disappeared and the 
ratio is the same for both distributions within 10 — 15%. 
The A dependence of the two parameterizations is quite 
different at lower Q^, however, quite similar at larger Q^. 

We now turn to the real part of the DVCS amplitude 
which depends both on the ERBL and the DGLAP re- 
gion. In Fig. 21 we plot the ratio of the real part of 
nuclear to nucleon DVCS amplitude in NLO for fixed 
and varying Xbj and in Fig.|31we show the same ratio for 
fixed Xbj and varying Q^. 

As can be readily seen in Fig.^jthe behavior of the real 
part, for both parameterizations, is substantially differ- 
ent from the one of the imaginary part for both parame- 
terizations except for the smallest values of Xbj ■ In partic- 
ular the dip in the ratio in the region 10"^ < Xbj < 10^^ 
is very surprising. The reason for this suppression is quite 
deep. To put it in a nutshell, it is the confluence of the 
nuclear medium effects: nuclear shadowing at small Xbj 
and anti-shadowing at large Xbj which mix under evolu- 
tion. How is that possible? 

The real part of the amplitude has two contributions: 
one originating in the ERBL region and one in the 
DGLAP region with a relative minus sign between them. 
If one were to decrease the contribution from the ERBL 
region and increase the contribution from the DGLAP 
region, one can obtain a significantly different value from 
the nucleon with just a relative small change in both con- 
tributions due to the relative minus sign. This is precisely 
what happens. To be definite let us take Xbj = 5 x 10~^. 
At this value of Xbj the value of the quark singlet is vir- 
tually not shadowed with gluon already experiencing a 
stronger suppression due to the relative shift to smaller 
values of X where nuclear shadowing is stronger. How- 
ever, due to the nuclear enhancement at larger X, the 
quark and gluon distributions in the ERBL are smaller 
compared to the nucleon case in order to fulfill the first 
and second moment of the nuclear GPD. At smaller Xbj 
nuclear shadowing leads to an additional suppression. All 
of this leads to a smaller value of the contribution to the 
real part from the ERBL region compared to the nu- 
cleon case. The contribution from the DGLAP region 
increases since the value of the contribution due to the 
nuclear enhancement increases. Therefore we have the 
required configuration where the ERBL contribution de- 
creases and the DGLAP one increases. This does not 
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FIG. 4: (Color online) Ratio of nuclear per nucleon to nucleon 
DVCS amplitudes vs. xi,j for two values of for the FGS 
(upper plot) and Eskola (lower plot) parameterizations. 
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FIG. 5: (Color online) Ratio of nuclear per nucleon to nucleon 
DVCS amplitudes vs. for four values of Xbj for the FGS 
(upper plot) and Eskola (lower plot) parameterizations. 



drastically change under evolution though the nuclear 
enhancement is pushed into the ERBL region the ac- 
tual evolution in the ERBL region is slower than in the 
DGLAP region. Furthermore, the nuclear shadowing in 
the giuon leads to a slower enhancement of the quark sin- 
glet compared to the nucleon case leading basically to a 
preservation of the effect for a larger range of except 
for the large values of Xbj ■ 

This can be observed in Fig. El where the ratio of the 
real part is plotted vs. for fixed xi,j and the suppres- 
sion effects is seen to decrease slowly at small Xbj and 
more rapidly at large Xbj though the absolute value of 
the effect remains large at large Xbj even for large Q^. 
Again, the slow change in the ratio is associated with the 
predominant enhancement of large light-like distances. 

This relative suppression of the ERBL region com- 
pared to the nucleon case suggests the simple picture 
that whereas in the DGLAP region partons could origi- 
nate in principle from different nuclei, in the ERBL this 
contribution is strongly suppressed. This once more has 
to do with the fact that the ERBL region at small Xbj is 
mainly sensitive to large light-like distances and correla- 
tions over large light-like distances are more susceptible 
to destructive nuclear medium effects leading to a sup- 
pression of these parton correlations. 

Since the suppression effect is very substantial, it 
should be clearly measurable in a high precision eA 
DVCS experiment at the EIC. We will discuss this in 
detail in a forthcoming paper. 

The above made observations would change somewhat 
if one were to consider non-zero t values. If we had per- 
formed a comparison at t = tmin = ~2;^jTO^/(1 — Xbj) 
rather than at < = the suppression of the nuclear ampli- 
tude at Xbj > 0.01 would have been even stronger (factor 
2 — 5 depending on the considered nucleus) since the nu- 
clear amplitude is suppressed by the nuclear form factor 
which can be approximated as Fa = exp(i?^t/6) for suf- 
ficiently small t, where Ra is the electromagnetic radius 
of the nucleus. 

In practical apphcations considered in the effects 
of the nuclear form factor for asymmetries are approxi- 
mately cancelled in the ratio. 

Furthermoew, for small Xbj nuclear shadowing effects 
are likely to lead to a small increase of the <-slope. How- 
ever, in the case of cross sections integrated over t (which 
is the only experimentally measurable quantity) this ef- 
fect leads to a reduction in the cross section by no more 
than 10% and it is even smaller for observables sensitive 
to interference terms. 

For Xbj > 0.1(200/^^/'^) deviations of the i-dependence 
from that of the nuclear form factor are much more 
sensitive to the details of the nuclear dynamics and re- 
quire a more microscopic treatment than the one pre- 
sented here. It should be remembered that for Pb — 208 
FA{tmin{x = 0.1)) is close to the point where the form 
factor changes sign and thus zero. Therefore, it will be 
extremely difficult to extract the coherent cross section 
for such an Xbj even using a veto for production of for- 
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ward neutrons. 

When comparing the two parameterizations one finds a 
substantial difference at low and again similar values 
at large except for the smallest values of Xbj where 
this difference persists over the entire interval. The 
most curious effect is seen at the smallest values of Xbj 
for the FGS parameterization where the nuclear medium 
effect is less for CA [Pb) compared to O (Pd). This 
is a peculiarity of the parameterization and once more 
the extreme sensitivity of DVCS to the gluon distribu- 
tion at small Xbj . Note that the Eskola parameterization 
does not have this feature since the gluon at small Xbj 
is substantially different from the FGS parameterization. 
Experiments will hopefully be able to distinguish this ef- 
fect in an observable sensitive to the real part like the 
beam charge asymmetry. 

IV. CONCLUSIONS 

In this paper we propose a model for nuclear GPDs 
based on the very successful GPD parameterization used 
in ^3 which describes all available world data of DVCS 



on the nucleon. We calculate the nuclear DVCS ampli- 
tudes using two different parameterizations for the nu- 
clear effects as proposed by Q and Q and plot the ratio 
of nuclear to nucleon DVCS amplitude for both the real 
and imaginary part. We find an enhancement of nuclear 
shadowing in the imaginary part for small Xbj compared 
to the inclusive case which can be traced back to the 
high sensitivity of DVCS to large light-like correlation 
distances which can be more easily influenced by nuclear 
medium effects. In the real part we observe for the first 
time an interplay between large and small Xbj nuclear 
effects in hadronic amplitudes for a large range in Q^. 
This interplay is most visible for 10~^ < Xbj < 10-1 
where the effect is strongest and should be observable in 
an eA DVCS experiment as for example at the EIC. 
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